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Abstract: We study the maximally supersymmetric BFSS model at finite temperature and 
its bosonic relative. For the bosonic model in p + 1 dimensions, we find that it effectively 
reduces to a system of gauged Gaussian matrix models. The effective model captures the 
low temperature regime of the model including one of its two phase transitions. The mass 
becomes for large p, with A the ’tHooft coupling. Simulations of the bosonic-BFSS 

model with p = 9 give m = (1.965 ± .007)which is also the mass gap of the Hamiltonian. 
We argue that there is no ‘sign’ problem in the maximally supersymmetric BFSS model and 
perform detailed simulations of several observables finding excellent agreement with AdS/CFT 
predictions when 1/a' corrections are included. 
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1 Introduction 

The leading candidate for a non-perturbative formulation of M-theory is expected to be the 
infinite matrix size limit of a matrix model of some kind. One such proposal is the BFSS 
model [1-3] which was conjectured to capture the entire dynamics of M-theory^. Relatives 
of this model such as the BMN model [4] or models derived from the ABJM model [5] are 
also considered possible viable candidates for such a non-perturbative formulation. All of 
these conjectured formulations of M-theory are regularised versions of the supermembrane 
and are matrix quantum mechanical systems. They are based on the matrix regularisation of 
membranes introduced by Hoppe [6] and extended to the supermembrane in [1] and [7]. 

In the second half of the paper we focus on the BFSS model. However, we dedicate the 
earlier sections of the paper to a careful study of the Hoppe regulated bosonic membrane 
which is the bosonic part of the BFSS model. This model is also of separate interest since it 
is the high temperature limit of maximally supersymmetric 1 -|- 1 dimensional SU (N) Yang- 
Mills (or matrix string theory) compactified on a circle when the fermions decouple due to 
their anti-periodic boundary conditions at finite temperature. The model has been studied 
already both theoretically [8] (see also ref. [9]) and numerically [10, 11] and our results are 

'^With a periodically identified lightlike direction for finite matrix size. 
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in accord with these studies. The zero temperature supersymmetric gauge theory has been 
studied in a hamiltonian light cone context in ref. [12]. 

We then continue our study with the BFSS model. This model first emerged as AA = 16 
supersymmetric Yang-Mills quantum mechanics [13-15], later it arose as the 11-dimensional 
supermembrane in Hoppe’s regularization and most recently as the BFSS model [2], a can¬ 
didate for a non-perturbative formulation of M-theory, and it also describes the low energy 
effective theory of a system of HO-branes [16]. A lattice version of the model which preserves 
eight of the sixteen supersymmetries was presented in [17]. 

Our lattice regularisation of the model follows Catterall and Wiseman [18]. In this 
formulation it is clear that the Pfaffian obtained when the Fermions are integrated out can 
in general be complex. However, we find that the phase of the Pfaffian is restricted to a 
narrow band so that cos(0pj) ~ 1. There is therefore no real phase or ‘sign problem’ as far 
as simulations of the model are concerned. We simulate the model using a rational hybrid 
monte carlo algorithm (RHMC) and find excellent agreement with results reported in [18-21] 
though our values for the energy are slightly higher than those of [19, 21] but in excellent 
agreement with predictions of AdS/CFT when leading Xjol corrections are included. Those 
interested primarily in the supersymmetric model can skip to section 4 for our discussion of 
the model and results, consulting section 2 for the continuum model and our notation. 

The principal results of this paper are: 

• We perform monte carlo simulations of the bosonic BFSS model and measure the two 
point correlation function, the mass gap and the eigenvalue distribution of each matrix. 
All fit with Gaussian matrix quantum mechanics with the same mass as that found 
from the correlation function. 

• We derive an effective description of the bosonic model using a 1/p expansion where p 
is the number of matrices. The description is in terms of p massive scalar fields that 
are gauged under the adjoint action of SU{N) but are otherwise free scalar fields^ 

• The effective model reproduces the one of two phase transitions (see [22]) of the full 
model with surprising precision^. 

• We study the maximally supersymmetric BFSS model and present arguments showing 
that though the lattice model in general has a complex phase, it is only the cosine of 
this phase that enters in simulations and the model is such that the angle is restricted 
to regions where the cosine is positive hence there is no sign problem in the full model. 

^The relevant results here were previously obtained in [8] which we became aware of these after writing this 
article. 

^The model can also be interpreted as the high temperature description of the maximally supersymmetric 
H-l dimensional Yang-Mills compactified on a circle. In this setting its two transitions are the high temperature 
limit of the black hole black string transition in the dual gravity model. 
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• We simulate the model using a Fourier accelerated rational hybrid monte carlo algo¬ 
rithm confirming results found earlier by other groups and find excellent agreement with 
predictions of AdS/CFT when subleading 1/a' corrections are included. 

The simulations we report provide a useful test of our code as we proceed to examine fur¬ 
ther observables and the inclusion of longitudinal M5-branes or equivalently Z14-branes. Our 
goal being the Berkooz-Douglas matrix model holographically dual to the backreacted D0/D4 
brane intersection. Such studies can provide a solid test of the AdS/CFT correspondence with 
flavour degrees of freedom, which is a widely used tool for non-perturbative analysis in flavour 
dynamics. 

The layout of the paper is as follows: In section 2 we introduce the BFSS model and 
continue in section 3 with a study of the bosonic part of the model describing our lattice 
discretisation and hybrid monte carlo algorithm. We then present our numerical results for 
this bosonic model and continue in section 3.5 to develop an expansion for the model in terms 
of the inverse of the number of matrices. This leads us to introduce our effective Gaussian 
model and compare it with the full model. In section 4 we return to the supersymmetric model 
and present our lattice discretisation of this model. We give a discussion of the Pfaffian phase 
of the model and present our results. The paper ends with a discussion of our results in 
section 5. 

2 The BFSS model 

The BFSS matrix model is a one dimensional supersymmetric Yang-Mills theory naturally 
arising in type IIA superstring theory as a low energy effective description of DO-branes. It 
is conjectured that in the limit of a large number of DO-branes, N, the model is equivalent to 
uncompactified eleven dimensional M-theory [2] while for finite N it corresponds to M-theory 
compactified on a light-like circle [3]. The easiest way to obtain the BFSS matrix model is 
via dimensional reduction of ten dimensional supersymmetric Yang-Mills theory down to one 
dimension. The resulting reduced ten dimensional action is given by [23]: 

Sm = ^ /^i^Tr{^(Do ^')2 + l[A^A^■] 2 -lT^Clo^ODo'I/ + iT^Clor[A^T]} , ( 2 . 1 ) 

where T is a thirty two component Majorana-Weyl spinor, F^ are ten dimensional gamma 
matrices and Cio is the charge conjugation matrix satisfying = —F^^. We take a 

representation for F^ in terms of nine dimensional (euclidian) gamma matrices 7 *: 

F* = 7 * (g) cTi , for i = I,..., 9 , 

F° = I 16 (g) ia2 , 

Cio = 6*9® ia2 , (2.2) 
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where Cg is the charge conjugation matrix in nine dimensions satisfying 6 * 97 *Cg = 7 * and 
(Tj are the Pauli matrices. With the following choice for the Majorana-Weyl spinor: 

^ = (0) ’ 

where il) is a sixteen component Spin{9) Majorana fermion and Wick rotating to Euclidean 
time, we obtain the Euclidean action: 

= , (2.4) 

which, as is manifest, is Spin{9) invariant. Note: we have not imposed any restriction on the 
nine dimensional spinor basis. Eor example if we choose 7 * to be in the Majorana representa¬ 
tion (where the nine 7 * are taken to be real and symmetric) then Cg = I 16 and we arrive at 
the more standard form for the action (2.4). However, we are interested in a basis in which 
the discrete theory is free of fermion doubling, which can be achieved by taking a basis [18] in 
which Cg = Ig ( 8 ) (Ti. We will return to this in section 4, where we consider the discretisation 
of the full BESS matrix model. 

3 Bosonic BFSS on the lattice 

In this section we focus on the bosonic part of the action (2.4) given by: 

Si, = . (3.1) 

The zero temperature model was introduced by Hoppe [ 6 ] as a gauge fixed and regulated 
description of membranes. The model has also been extensively studied in the literature^. It 
has been simulated for a first time in refs. [24, 25], where certain aspects of the full model 
were described in terms of simple Gaussian model. Its phase structure at finite temperature 
has been explored in [ 8 , 10 , 22 ], where the authors found that as the temperature is decreased 
the model first undergoes a 2 nd order deconfining-confining phase transition into a phase 
with non-uniform but gapless distribution for the holonomy. As the temperature is further 
decreased there is a 3rd order transition to a gapped holonomy with a quadratic decrease in the 
internal energy to a constant value for lower temperatures. The high temperature expansion 
of the model was developed in [26]. In what follows we will reproduce the main results 
of [ 8 , 22] and elaborate on the properties of the theory at zero temperature. In particular we 
will provide evidence that the low temperature phase of the model has an effective description 
in terms of a free massive scalar which captures many of the finite temperature features of 
the model including one of its two phase transitions. 

■^The model is also the high temperature limit of 1 -I- 1 dimensional Af = 8 supersymmetric Yang-Mills on 
R X where jS is now the period of the 5^ and the fermions drop out due to their anti-periodic boundary 
conditions at finite temperature. 
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3.1 Discretisation 

We discretise time to A sites tn = an, (n = 0 ,..., A — 1), where the lattice spacing is a = /3/A, 
and impose periodic boundary conditions identifying the point t\ = Aa = /3 with the point 
0. To discretise the covariant derivative 'Dt we define the transporter helds: 


Un,n+l = v exp 


'*(n+l)a 


dt A{t) 


(3.2) 


where V denotes a path ordered product. Let us consider for a moment the pure derivative 
part of Df On the lattice we have: 


dfXl 


K+i 


- XI 


(3.3) 


To make the above expression gauge covariant we have to transport back the field at to 
tn- For the discrete version of the covariant derivative, we obtain: 


Vt 


Un,n+lX'!^j^lUn+l,n — X\, 


(3.4) 


where Un+i,n = Ulm+i- Using equation (3.4) for the discrete bosonic action we obtain: 


A-l 


Sb = N'^ti 


n=0 


— Ku. 

a 


n,n+lXl^lUl 


+ - {X^f 



(3.5) 


where without loss of generality we have taken g = The action Si, can be written in a 

much simpler form by using the U(n)^ gauge symmetry of the model. Indeed, at each lattice 
site we have a local U(A^) symmetry. Using that symmetry we can perform the transformation: 


K = W* , (3.6) 

X'\ = Uo,iXiul, , 

X'i = {Uo,iUi,2) Xi{Uo,iUi,2y , 


^1-1 = (Uo,iUi,2 ... Ua-2,a-i) X{_, (C/o,iUi,2 ... Ua-2,a-i)^ 

introducing the notation W = (?7o,if7i^2 • • ■ Ua- 2 ,A-iUa-i,o) for the bosonic action (3.5) we 
obtain: 


= _liVtr I ^ X'lXl^, + X'\_,W X'lW^ I + AT ^ tr I - (W^)^ - ^ [X'l X^f] 
® U=0 J n=0 J 


(3.7) 


“This can always be arranged by an appropriate rescaling of the matrices and the time coordinate and fi 

, 1/3 2 

becomes the dimensionless temperature parameter with X = g N the’t Hooft coupling. 


- 5 - 







The unitary matrix W has the decomposition W = VDV\ where D = diag{ • • ■) is 
a diagonal unitary matrix and 1/ is a unitary. But the action (3.7) has the residual symmetry 
X'l VX'lV\ which we can use to diagonalise W. Furthermore, it has an additional 
symmetry X'\ —)• hnX'\^hh, where is a diagonal unitary matrix, which we can use to 
“distribute” the diagonal matrix D among all of the hop terms. Indeed, defining the matrix 
Da = diag{e®®i/^,..., which satishes {D\)^ = D, one can verify that under the 

transformation: 

X'l = {VK)K{VK)^ , where : K = (Da)- , (3.8) 

the action (3.7) transforms into: 


A-i 


5fe[X,DA] = J^tr 


n=0 




i \2 



(3.9) 


Now let us discuss the measure of the transporter helds The measure can be written 

as: 

A-l A-l A-l 

nw„ ,n+l — n DD „,„+1 DDo,i = n VUn,n+i vw , (3.10) 

n=0 n=l n=l 

where we have used that Do,i = W (Di ,2 • ■ ■ Da- 2 ,A-i)^ and the invariance of the measure. 
But the action (3.9) depends only on the matrix W (in fact only on the eigenvalues of W). 
Therefore the integration over the measure of the transporter fields reduces to: 


. A-l 


oc (3.11) 


n=0 

. N . r ^ 

oc / n ^-Sblx,D/,{e)] ^ 

k=l l>m k=l 


YldOk 


,-Si,[X,DA(6l)]-SFp[0] 


where S'fp[^] is given by: 


S'fp [ 6 ] 


= -E>» 

l^m 


. 0i — 9m 

sm- 

2 


(3.12) 


3.2 Hybrid Monte Carlo 

To implement the Hybrid Monte Carlo algorithm it is convenient to work in a gauge in which 
the holonomy matrix is non-trivial at only one link (we choose the one connecting sites zero 
and A — 1). To this end we omit the diagonal matrices hn in the transformation (3.8). The 
action (3.9) is then given by: 


56[X, D] = iVtr <^ — 


f A-2 A-l 

E + E 

I n=0 n=0 


i \2 


-{K) 

a 
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The corresponding Hamiltonian for the molecular dynamics part of the HMC algorithm is 
then: 


nb 


os — 



(3.14) 


where are the canonical momenta corresponding to the hermitian matrices and are 
the canonical momenta associated to the angles 9i. Hamilton’s equations read: 


Kin. = -dSb/dXl^^i , = -dSb/dOi , (3.15) 

vi = P* R, = 

^n,lm ^n,lm ) ^l i 


where dots denote derivatives are with respect to the Monte Carlo time. Using equation 
(3.13), for the corresponding forces we obtain: 


N 


a 

N 


-dSb/dXl^i = -(x{-2Xl + D^Xl_M^ +Na 


Im 




xi.xi 


- J Im 


-dSb/dXl^i = -{xi^,-2Xi, + Xl_,)^^ + Na[Xi,[XlXi]]^^ for n = l,...,A-2, 
/ ^ ^ \ 

^A-V ^A-l^^A-1 


-dSb/dXX.,,mi = - [DX^oD^ - 2Xl_i + + Na 

d \ / Im 


- J Im 


2N 


N-l 


m=0 




m , mj^l 


Ol-Or, 


(3.16) 


which we implement in the Hybrid Monte Carlo. 


3.3 Phase structure 

In this section we reproduce the studies of the phase structure of the bosonic model originally 
obtained in [22] . The main observables that we focns on are the internal energy of the system 
E, the “extent of space” (P^) and the Polyakov loop P dehned via: 


/ ^ \ 

( iNp J ■ 

(3.17) 

/ 1 r \ 


' 0 ' 

(3.18) 

“a 

III 

(3.19) 

/ ? \ 


U = V exp i i dt Ao{t) \ , 

(3.20) 


0 


where the holonomy matrix, U, is the continuum limit of the link variable ?7o,A defined in 
equation (3.2). The expectation value of the Polyakov loop (jPj) plays the role of an order 
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Figure 1. A Plot of the expectation value of the Polyakov loop |P| as a function of the temperature, 
for N = 16 and lattice spacing a « 0.05. One can see that near T « 0.90 the theory undergoes a 
second order phase transition. 

parameter for the confining-deconfining phase transition discussed in [22]. In figure 1 we 
presented a plot of this order parameter as a function of the temperature. The plot is for 
= 16 and lattice spacing a ~ 0.05. One can see that near temperature T ~ 0.90 there is 
a second order phase transition. The change of the slope of the curves and the fluctuations 
in the simulations near T ~ 0.9 is consistent with the existence of a second order phase 
transition. 

In figure 2 we present plots of the energy and “extent of space” as functions of the 
temperature, for = 16 and lattice spacing a ~ 0.05. The dashed curves represent the 
analytical high temperature behaviour obtained in [26] . Our results agree very well with the 
corresponding studies in [8, 10, 22]. 

A more detailed analysis of the temperature range close to the transition revealed that 
there are in fact two transitions. To uncover more detail on the nature of the phase transition 
the authors of [22] analysed the distribution of the holonomy matrix near the phase transition 
and uncovered behaviour consistent with the Gross-Witten model [27] and concluded that the 
holonomy undergoes a transition from a uniform distribution at Tc 2 — 0.8758(9) to a gapped 
distribution at Tc 2 = 0.905(2). In figure 3 we present our plots of the distribution of the 
holonomy around the phase transition. The dashed curves in the plot represent fits with the 
gapped and ungapped forms of the Gross-Witten distribution which are in excellent agreement 
with those of [8, 10, 22]. We have not attempted to refine their results, rather our purpose is 
to uncover a hidden Gaussian structure in the model. 
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Figure 2. Plots of the scaled energy EjN"^ and the “extent of space” (i?^) as functions of the 
temperature. The dashed curves correspond to the high temperature behaviour obtained in [26]. One 
can see that near T « 0.9 the plots suggest the existence of a second order phase transition. The 
energy and temperature in the plots are in units of 



e 


Figure 3. Plots of the distribution of the holonomy P for temperatures T = 0.900 (the confined or 
ungapped) and T — .9006 (the deconfinied or gapped). The plots are for size iV = 16 and lattice 
spacing a « 0.05. The dashed curves correspond to fits to the Gross-Witten gapped and ungapped 
distributions. 


3.4 Gap and eigenvalue distribution 

In this section we investigate the eigenvalue distribution of the scalar fields. We also perform 
studies of the mass of the theory at zero temperature. Our results suggest that at all tem¬ 
peratures the eigenvalue distribution of any one of the is given by a Wigner semicircle, 
with a radius R\ following the temperature behaviour of the observable Therefore, we 

conclude that while the radius of the distribution experiences a phase transition the shape of 
the distribution remains unchanged. We believe that the main reason for this behaviour is 

®The semicircle law implies 
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that for nine scalar fields the commutator squared term can be replaced and approximated 
by a quadratic mass term in the spirit of [28], where an expansion at large number of scalar 
fields has been developed. Generalising these techniques, we are able to obtain an estimate 
of the mass, which agrees very well with both the gap measured from correlation functions 
and the radius of the distribution which are obtained from Monte Carlo simulations. 

In the limit of high temperature the model reduces to the 10-dimensional Yang-Mills 
matrix model considered in [28]. The obtained behaviour of the radius is in agreement with 
the large temperature expansion performed in [26] and provides an analytic approximation 
to the leading coefficients in this expansion. 

We begin by considering the model at zero temperature. In this case the holonomy can 
be completely gauged away and the model simplifies. Furthermore, at zero temperature the 
correlator: 

(Tr(Y^(O)Y^(t))) (xe-™‘ + ... , (3.21) 

captures the gap Am = Ei — Eq of the theory. To calculate the gap in the discrete theory, 
we periodically identify the time direction with period /?: 

{Tv {X\0)X\t))) = A{e-^^ + , (3.22) 

Note that although formally /3 is the same parameter that we have at finite temperature, 
since we set the holonomy to zero here its meaning is just a periodic coordinate as opposed 
to inverse temperature. Our result for the correlator for N = 30, /3 = 10 and lattice spacing 
a = 0.25 is presented on the left in hgure 4. The fitting curve is given by equation (3.22) and 
when we perform a two parameter fit we obtain A ~ 7.50 ± 0.2 and m ~ (1.90 ± .01) . 

However, for Gaussian scalar fields of mass m we have A = performing a 

one parameter fit for m yields m = 1.965 ± 0.007 and A = 7.63 ± 0.03 On the right we have 
presented a plot of the eigenvalue distribution of one of the matrices for the same parameters. 
The fitting curve represents a Wigner semicircle of radius Rx ~ 1.01. The fact that the 
theory is gapped and that the eigenvalue distribution is a semicircle suggests that that at low 
temperate the model has an effective action: 

OO 

S’eff = A^ I dtTr Q(W)2 + lm2(W)2^ (3.23) 


for each of the matrices X*. It is well known [29] that for the action (3.23) the eigenvalue 
distribution of X is given by a Wigner semicircle of radius: 





1.009 ±.002, 


(3.24) 


where we have substituted m ~ 1.965±.007. This agrees nicely (within errors) with the result 
for Rx ~ 1.01 obtained by htting the actual distribution. It is also in excellent agreement 
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Figure 4. On the left: A plot of the correlator (Xr for N = 30, /3 = 10 and lattice 

spacing a = 0.25. The fitting curve is given by equation (3.22) with A = with 

parameters m « (1.965 ± .007)A^/^ . On the right: A plot of the eigenvalue distribution of one of the 
scalars for the same parameters. The fitting curve represent a Wigner-semicircle of radius R\ « 1.01. 


with the large p theoretical prediction of [8], 

+ ^(^ - |) + ■ ■ ■) - 1.0068 (3.25) 

3.5 1/D expansion 

Adapting the techniques developed in [28] to the time dependent case that we are considering 
we can obtain a theoretical estimate of the radius at low temperature^. Let us consider again 
the action of the bosonic model (3.1): 

Sb = N , (3.26) 

where we have rescaled so that ft = (effectively we set g = The commutator 

square term can be written as: 

Tr[X\ = Tt ([t“, ft] [t\ t'"]) Xixl X/Xi = X/X^ , (3.27) 

where D are SU{N) generators normalised to TiDt^ = <5“^ and the tensor is given by: 

It is convenient also to define the inverse kernel of satisfying: 

Rabcd = 6l6l , = 5/5) . (3.29) 

^The next order corrections for the current model were computed in [8] but we were unaware of this work 
at the time of writing. 
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We will also use the identities: 


X-bcd ^_2N , ^iabcd = ^ab • 


(3.30) 


The action (3.26) can then be written as: 


Sb = N dtTi ( ^(AW)^ ) - ^—X 


1 \ AT 

^^2'\ , abed I dtXlXlXiXy 


10 


Our next step is to add to the action the term®: 


A5 = J^^abcd I dt + X'^^^’^XlxCj , (3.32) 


/3 


(3.31) 


the action Si = Sb + AS then becomes: 


„ N 




^ dt {tV(AW)2 + k^'^xixi'j + Jtiabcd I dtkV . (3.33) 


kiUr.^ = tf,k<^HL , (3.34) 


Next we define 

’ ij,lm — ^ij T^lrn ) 

From the definition of kij^im it follows that it is traceless with respect to the first and the 
second pair of indices and we can invert: kji^mi Substituting in the action (3.33), 

Fourier transforming (via X = e Xn) and assuming kab is time independent we obtain: 


N 


= v E 


2 ^ \ V 0 

n=—oo \ ^ ^ 


27rn + 0,-0A _ , , _ \ ~ , /3iV _ 


Pjl,ml + kjl,ml]Xn,lm+ ^ ^abedk^^k^ , (3.35) 


where we have defined the projector on traceless matrices: 

1 


Pij,lm ^ij ^Im dirndjl j^dij6ijn 
and assumed that k is constant. Defining also the double index matrix : 


IF {jl)ij^lm — 

we can integrate out the X’s: 


27rn + — 6j 

~0 


Pjl,ml T kjl,ml ; 


j VXVke-^'b (X j Vk JjDef^/^ {P.W{n).P) 


g-^b^abcdk^'’k‘= 


(3.36) 


(3.37) 


(3.38) 


^Note that we can always add this term since f Dke ^ = const 
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The effective action for the held k then becomes: 


^efr[A:] = f E Tr log {P.W{n).P) + ^^^abcd kV 


(3.39) 


We now notice that the hrst term in the expression for the matrix W (3.37) commutes with 
the projector P. It is natural then to consider an ansatz for k which also have that property. 
Thus we consider: 

kij^lm — kij Pij^lm — Pij,lm klm (3.40) 

The last equality is possible only if all diagonal components of kij are the same (namely 
kii = kjj for all i,j) we also choose kij to be symmetric. Then: 


Pij,lnn Pij,lm^lmip) t 


(3.41) 




27rn + 9i — 9j 


T kij , 


and we have for all powers r that {P-W{n).P)^j = Pijjm ■ Therefore, 

Ses[k] = Tijji log(A(?T,)jj) + l^abcdk k 




The corresponding saddle point equation = 0 becomes: 


E 


p 




+A:, 


I ^ I„a,b+C +d _ n 

T 2 l-^abcd ^ ^ 


(3.42) 


(3.43) 


We can now sum the series to obtain: 

p sinh(/3Y^) 


‘2‘NIdabcdk^'^tijt'j^ — 0 


(3.44) 


cosh{/3^/kY) - cos(9i - 6j) 

In principle we can solve for kij in terms of 9i — 9j , but we will restrict ourselves to extract the 
low temperature dependence. The hrst term in equation (3.44) has the following expansion: 

sinh(/3i/%) 


cosh.{l3 ^Jkij) — cos{9i — 9j) 


= 1 + 2cos{9i - 9j)e + 0(e 2/3./^) 


(3.45) 


One can see that the effect of the holonomy is exponentially suppressed at low temperature, 
To leading order we can then consider a symmetric ansatz kij = ko, which also implies 
= ko We obtain: 


-^P 


ko Pij,ji — 0 • • ko — p 


— 


Substituting into the action (3.33) to leading order we obtain: 

^ / 2/3 ^ 

S, = N f diTr f i(X')2 + ^(x')2 


(3.46) 


(3.47) 
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Figure 5. On the left: A plot of the energy of the gauged gaussian model for N = 32 lattice spacing 
a « 0.05. The red curve represents Wigner semicircle. On the right: A plot of the expectation value 
of the Polyakov loop for N = 32 lattice spacing a « 0.05. The second order phase transition takes 
place at T « 0.95. The energy and temperature in the plots are in units of 


and the corresponding radius of the eigenvalue distribution is: 

Rq = 1 - j 0.98 . (3.48) 

This result agrees within a few percent with Rx ~ 1.01 obtained by fitting the distribution 
in figure 4. 

The exponential suppression of the holonomy corrections to the low temperature saddle 
point for suggest that a lot of the physics of the model (at least at low temperature) can 
be captured by the action: 


Sb = N 


j dtTr ( ^(AX')2 + 
0 




(3.49) 


where we restored the gauge field in the covariant derivatives. Surprisingly this model knows 
lots about the phase transition of the full model. An analysis similar to that above shows that 
the critical temperature for p adjoint gauged Gaussian matrices as in [9] with mass m occurs 
at RGaussian _ ^ which for m = yields Tc = 0.9467. In figure 5 we have presented 
our results for the energy and the Polyakov loop. Note that in this approximation the scaled 
energy E/N"^ is equal to the extent of space {E?). The plots are for N = 32 and lattice 
spacing a ~ 0.05. One can see that both the energy and the Polyakov loop exhibit the same 
behaviour as for the full model. There again appear to be two distinct transitions with the 
higher temperature one taking place at T ~ .95 and again appearing to be second order. It 
is slightly shifted towards high temperatures relative to the critical temperature for the full 
model. If instead of mass m = we use the value m ~ (1.90 ± .01)A^/^ the phase 

transition is shifted in the opposite direction (just bellow T = 0.9). This indicates that if one 
fits the mass parameter one can improve even further the agreement of the gauged gaussian 
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Figure 6. A zoomed in plot of the energy versus the temperature. The dashed curve represents a fit 
to equation (3.50) with Tg « 0.896. 

model and the full one. The dashed curve in the second plot is the theoretical prediction 
for the high-temperature behaviour of the Polyakov loop, again one can observe a very good 
agreement. The high temperature behaviour of the energy on the other side disagrees with the 
corresponding behaviour of the full model. This is not surprising since we derived the effective 
action at low temperature and the dominant behaviour at high temperature is dominated by 
the highest power of the potential which has been changed from quartic to quadratic. One 
can also see that at low temperature the energy remains constant. 

In figure 6 we have presented a plot of the energy versus the temperature zoomed in near 
the phase transition. The dashed curve represents a fit with: 

E-eo = C{T- Tof (3.50) 

with parameter Tq ~ 0.896. This indicates that the third order phase transition that the full 
model exhibits [8, 22] is also captured by the gauged gaussian model. 

One can perform the large p analysis (see [28]) in the high temperature limit where the 
model now has p+1 matrices (the holonomy becomes the additional matrix) and predict that 
the model in this limit again becomes Gaussian but now the field kab includes the holonomy 
and the saddle becomes kab = Vp + ^^ab which again predicts a Gaussian matrix model 
with a high temperature effective action Sgff = and consequently a Wigner 

semicircle distribution for the eigenvalues os X^. 

We conclude this section by presenting results for the eigenvalue distribution of the gauged 
gaussian model at hnite temperature. In figure 7 we presented plots of the distribution for 
T = 0.2, N = 30 (left) and T = 5.0, N = 30 (right). The red dashed curves show a Wigner 
semicircle. One can see that the shape of the eigenvalue distribution does not change with 
temperature. 
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Figure 7. On the left: A plot of the eigenvalue distribution for N = 30, T = .2 and lattice spacing 
a = .25. The red curve represents Wigner semicircle. On the right: A plot of the eigenvalue distribution 
for N = 30, T = 5.0 and lattice spacing a = .05. Again the red curve represents Wigner semicircle. 


4 The supersymmetric model on the lattice 

In this section we consider the full supersymmetric BFSS model on the lattice. The model has 
been simulated using non-lattice approach in [19] and using lattice discretisation in [18] and 
[21]. The main goal of these studies has been to compare the low temperature regime of the 
model to the holographically dual black hole geometry. The authors of refs. [19] and [21] also 
compared the high temperature regime of the model with the high temperature expansion 
performed in [26]. Our goal is to reproduce some of these studies and calibrate our code. 

A naive discretisation of the action (2.4) would result in a fermion doubling. This can 
be avoided [18] if the charge conjugation matrix is taken to be Cg = Ig < 8 ) cti. ® Constructing 
a basis for which Cg is of this form is relatively straightforward. For example one can tensor 
up the Majorana basis in seven euclidean dimensions 7 ^: 

7 “ = - 7 I < 8 > 0-3 , for a = 1,..., 7 , 

7 ® = Is < 8 ) fT 2 , 

7® = Is (8) (Ti , (4.1) 

and verify that indeed Cg is of the desired form (it also satisfies Cg = 7 ®). We next proceed 
in discretising the action (2.4). Since the bosonic part of the action is identical to the one 
considered in section 3 , we will consider only the fermionic part of the action: 

5/ = ^/drtr{v^“Cg„^P,V’^ -V’“(Cgy)„^[X',V'^]} . (4.2) 

We begin by splitting the fermions into two eight component fermions: ^|J = (V'i,V' 2 ) and 
®It is analogous to using staggered fermions, which in one dimension complete removes the doublers. 
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defining the forward and backward derivatives D±\ 


{V_W)n = [Wn - Un,n-lWn-lUn-l,n)la , 

{V+W)n = {Un,n+lWn+lUn+l,n " Wn)/a . (4.3) 


One can show that the discretised kinetic term then becomes: 

1 r 

•S'/” = —^ dr tr ^ ^ tr {V’?^„(P_V’ 2 )n + V’i^n(^+V'l)n} = (4.4) 

^ ^ n=0 

^ r A-1 A-2 

= \ ~ '4’2,nUn,n+lill,n+lUn+l,n ± V’2iA-lt^A-l,0^1,0^0,A-l 

^ I n=0 n=0 

where the plus/minus sign in the last term corresponds to periodic/anti-periodic boundary 
conditions for the fermions.^® Using the gauge from the previous subsection when the holon- 
omy is concentrated on a singe link we can write S'j™ as: 



ckin 




{ A-l A-2 

n=0 n=0 


n+l ± V'i’A-l-^^V’1,0 D"' 


(4.5) 


Since all fields transform in the adjoint of SU{N) instead of dealing with matrices we can 
use the corresponding real components: = tr(t“X) and = tr(t“?/;), where are the 

standard basis of SU{N) (introduced earlier) normalised as trt“t^ = 5“^. Sj™ can then be 
written as: 


Af2-1 A-l 8 

■S/" = 3 E E ■ ('‘^<5) 

a,6=0 m,n=0 q;=1 

^mn = (<5m+l,n “ <^m,n)<J“^ ± <5m,A-l<5n,0 (4.7) 

. (4.8) 


where the plus/minus sign corresponds to periodic/ant-periodic boundary conditions. The 
kinetic term can also be written as: 


ckin 


M 


kin 

mn^a0,ab 


N^-l A-l 16 


a,6=0 m,n=0 a,/5=l 

Os -Ktr 
2ff2 Os 


OLfB 


Discretising the potential part of the action is straightforward. One obtains: 


(4.9) 

(4.10) 


N^-l A-l 16 




«/“= E E E 

a,6=0 m,n=0 a,/3=l 


(4.11) 

(4.12) 


^^Naniely the conditions ip-i = ±ipA-i and ipA = dzipo. 
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Finally, defining: 


M 


mn,a/3,ab 


Os -Kt, 

2g^ [k± 08 


+ {CgJ%pX^/6n+mfl 

ap 


(4.13) 


We can write: 

Sf = ifP"Mij) . 


(4.14) 


4.1 The PfafRan phase is not a problem! 

Integrating out the Fermions, the partition function of the model can be written as: 

Zoc j VX V'lp (X ypXPf(AI)e"'^‘>“[^l (4.15) 

Observe that A4 is the sum of an anti-hermitian kinetic term and a hermitian potential and 
Xi^{X) = —Xi{—X). Also because the Bosonic action is symmetric in X the PfafRan in the 
partition function can be replaced by |Pf(AI)| cos(0p/). Now as long as —^ < 0p/ < | the 
cosine is positive and the effective action defines a true probability distribution given by 

Seff = 5bos[W] -ln|Pf(A4)| -lncos(0p/) (4.16) 

In figure 8 we presented plot of the phase of the pfaffian of the fermionic matrix for N = 3 
and four lattice points.One can see that the cosine remains positive. We believe that the 
drop in the cos0 curve at very low temperatures is a lattice effect and is not present in the 
continuum limit. Our results show a very good agreement with the earlier studies of ref. [30]. 

4.2 RHMC and fermionic forces 

The next step is to apply the RHMC method [31] to the model. To this end we need the so 
called fermionic forces. Let us summarise briefly the philosophy. 

As we have shown above the model does not suffer from a severe sign problem and so we 
ignore the phase of the PfafRan and use that: 

lPf(M)l = det(Mt, (4.17) 

to write 

Zoc / , ( 4 .I 8 ) 

where 

Vf = (4-19) 

^^Note that to control the flat directions at low temperature we have added a small mass for the bosonic 
held. 
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Figure 8. A plot of the phase of the pfaffian of the fermionic matrix for N = 3 and A = 4. One 
can see that the phase remains small for all temperature, but drops at ver low temperatures. We 
believe that this is a lattice artifact and is not present in the continuum limit. The flat directions are 
controlled with a small mass regulator at low temperature. 



Here ^ is a 16(iV^ — 1)A dimensional vector consisting of the pseudo-fermionic fields. The idea 
of the RHMC is to approximate the rational exponent of the matrix Ai with a partial 
sum: 

# 

= ao + Y,<^dM^M+Pi)-\ (4.20) 

i=l 

where the parameters ao, Oi, (3i and # depend on the rational exponent d, the spectral range of 
the matrix Ai) At and the desired accuracy. We will need two rational exponents. To update 
the pseudo fermions we use that the held rj = {At^ At)~^/^^ has a gaussian distribution 
and solve for ^ = {At^ At^^^ rj using a multi-shift solver. Therefore, <5 = 1/8 is one of the 
rational exponents that we need. To calculate the fermionic forces and the contribution to 
the hamiltonian we need to invert (Al^ At)~^^^ and the second exponent is (5 = —1/4. 

Let us elaborate on the computation of the fermionic forces. We have two type of forces: 
gradients with respect to Xnij and gradients with respect to the phases of the links 9i. 
Using the partial expansion (4.20), one can easily derive expression for the derivatives of the 
fermionic action: 


<9Sps.f 

du 


- M + A)-' {M^ M + A)-'C 


# 

h. 

i=l 


t d{A[^ At) 
du 


(4.21) 


where hi satisfy (Al^ At + /?,) hi = ^i and are obtained from the multi-solver. 
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4.3 Simulation results 


In this subsection we provide our results from the Monte Carlo simulation of the model. We 
focus on the same observables that we analysed for the bosonic model in section 3.3. The 
definitions of the extent of space and the expectation value of the Polyakov loop P 

remain the same. The expression for the internal energy is [18]^^: 

= ~Jp (^('S'bos) - 2 ^ “ 1)^ > {4:.22) 

We have simulated the following configurations. For temperatures T > 2 we have used N = 8 
and A = 8 . For the region 1 < T < 2 we have used N = 8 and two or three different sizes 
of the lattice (for each point) in the range 8 < A < 16 (For T = 1 we also went to A = 32). 
For temperatures lower than one we have used N = 10,12,14 and two lattice sizes per point 
A = 8,16. For all temperatures the Polyakov loop is largely independent on the lattice 
spacing. The extent of space also experience very weak lattice effects. However, this is not 
the case for the internal energy and even for temperatures as high as T = 2 lattice effects can 
be a factor. In figure 9 w present our results for the energy at T = .9,1.0 for different lattice 
spacing. One can see that the lattice effects die out linearly, which allows us to extrapolate 
the energy to zero lattice spacing. Our results for the internal energy are summarised in 
figure 10. The dashed curve at high temperatures is the theoretical curve obtained in the 
high temperature expansion. The dashed curve at low temperature represents the prediction 
of AdS/CFT.^^ The model becomes unstable for small matrix sizes N, an effect which has 
been related to Hawking radiation in the dual gravitational theory [30, 32]. To compare 
with the AdS/CFT predictions one needs to consider large matrices. Simulations with large 
matrix sizes is computationally expensive and as a result our low temperature data is still 
preliminary. However, even at this point we have excellent overall agreement with the studies 
of [19] and [21]. 

Finally, in figure 11 we present our results for the Polyakov loop and the extent of space. 
Again the dashed curve represents the high temperature theoretical result obtained in [26]. 
One can see the excellent agreement at high temperatures. Our result for these observables 
agree with the previous studies preformed in [19], [18], [20] and [21]. 

5 Discussion 

In this paper we have analysed both the purely Bosonic and the supersymmetric BFSS models. 
These are “Hoppe” regulated membranes and one would expect that in the large N limit when 

^^Note that this expression is also valid for the bosonic model. The result can be obtained by rescaling the 
fields such that the kinetic term is temperature independent and removes any temperature dependence from 
the measnre (the Van Vleck Morette determinant generically depends on temperatnre). Then differentiating 
with respect to temperature and using the Ward identities associated with the total number of degrees of 
freedom yields this result. 

^^We have used the a' corrected expression = 7.41 ( j — 5.58 ( ) '' obtained in [20]. 
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Figure 9. Scaling of the internal energy with the lattice spacing for temperatures T = 0.9,1.0 and 
N — 10. One can see that the lattice effects dies our linearly, which allows extrapolation of the zero 
lattice spacing result. 
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Figure 10. Results for the internal energy for 8 < iV < 14 and 8 < A < 16. The dashed curve at high 
temperature correspond to the theoretical results of [26], while the low temperature curve represent 
the prediction for the internal energy from the AdS/CFT correspondence. 


the regulator is removed that they describe the full quantum dynamics of these membranes. 
Surprisingly we found that the bosonic model, for sufficiently large embedding dimension 
reduces to a system of p-massive free bosons with the mass given by m ~ . For p = 9 

we performed detailed simulations of the model evaluating both the fall off of the correlation 
function and the eigenvalue distribution of the X* fit with a Wigner semi-circle both of which 
give a consistent mass m ~ (1.965 ± .007)A^/^. This is a fundamental non-perturbative result 
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Figure 11. Plots of the expectation value of the Polyakov loop (|P’|) and the extent of space (i?2) 
as functions of temperature. The dashed curves represent the predictions of the high temperature 
expansion. 

and gives the mass gap in the full Hamiltonian of the model. 

The correspondence of the full and gauged Gaussian model is excellent for a wide range 
of temperatures. Somewhat surprisingly the phase transition region of the full model is 
faithfully reproduced by the effective model with the two transitions of the full model merged 
into one. Since the finite temperature version of the model is also the high temperature limit 
of 1 + 1 dimensional maximally supersymmetric Yang-Mills [24] compactified on a circle, we 
have established that this latter model should also reduce to a system of free massive scalars 
in its large radius high temperature phase. 

We then study the full supersymmetric BFSS model using a rational hybrid monte carlo 
simulation with Fourier acceleration to evaluate observables of the model. After describing 
our lattice discretisation of the model we investigated the phase of the Pfaffian obtained on 
integrating out the Fermions. The Pfaffian is generically complex, however, its phase is in 
fact not a problem for simulations. What enters simulations is the cosine of this phase and 
in the regularisation used in our work this phase is in fact restricted to a region where the 
cosine is positive once the lattice spacing is sufficiently small. Direct measurements confirm 
that the phase is indeed small. 

Though our results for this part of the paper do not yet go beyond those of [19] or cover as 
low a temperature as those of [21] they are more precise than those of Catterall and Wiseman 
[18] who use a similar lattice simulation. We have taken several lattice spacings and then 
performed an extrapolation to the limit of zero lattice spacing. We find good agreement 
with earlier results and excellent agreement with the predictions of AdS/CFT once Xjol 
corrections are included. Our results appear to approach the predictions of AdS/CFT a 
little more closely than those of [19] but the difference is broadly within the errors. The 
principal difficulty in simulating the model at lower temperatures is due to critical slowing 
down and though Fourier acceleration helped for T ~ 0.5, simulations becomes increasingly 
more difficult at lower temperatures. An additional difficulty is the instability due to flat 
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directions, which require increasingly larger matrix sizes as the temperature is reduced. 

One of the principal aims of this work was to check the claims of previous work and 
in particular those on the absence of a complex phase problem. We were also interested in 
calibrating our code as we extend it to include systems with L)4-branes. The extension to 
such systems will allow us to perform more extensive tests of the AdS/CFT correspondence. 

Acknowledgements: We wish to thank S. Catterall, M. Panero and M. Hanada for 
helpful comments. 
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